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Abstract 

We theoretically investigate quantum fluctuation of charge between even and odd states of a normal- 
superconducting-normal single-electron tunneling transistor. It is shown that due to the superconducting gap, the 
charge fluctuation in the Coulomb blockade regime for even state is larger than that for odd state. We show that 
large energy correction in the former regime caused by charge fluctuation can be explained by considering the 
charging energy renormalization. 
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The charge fluctuation caused by strong tun- 
neling is one of the basic problems in physics of 
single-electron tunneling (SET) transistors. Re- 
cently, there has been much development in the 
theoretical research [1-3]. Schoeller et al. calcu- 
lated the current by using a systematic diagram- 
matic technique based on the Keldysh formalism 
and predicted the renormalization of conductance 
and charging energy. The conductance and the 
charging energy decrease with decreasing tem- 
perature or bias voltage. Recently, Joyez et al. 
observed the conductance renormalization at low 
temperature [4]. 

As for superconducting islands, the parity ef- 
fect[5,6] has been observed. The parity effect is an 
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evidence that a Cooper pair is composed of two 
electrons. Though there is an investigation on the 
large charge fluctuation which takes account of the 
parity effect [7], there is no study which focuses on 
the quantum fluctuation of charge between even 
and odd states. 

In this paper, we mainly study the ground-state 
properties of a normal-supcrconducting-normal 
(NSN) SET transistor with a special attention 
to the charge fluctuation between even and odd 
states. This paper is complementary to a previous 
paper [8] where we studied transport properties. 
We assume that the charging energy Ec exceeds 
the superconducting gap A. We show that the 
charge fluctuation in the Coulomb blockade (CB) 
regime for even state is larger than that for odd 
state. The large energy correction caused by large 
charge fluctuation in CB regime for even state can 
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be explained by considering the charging energy 
renormalization. 
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Fig. 1. (a) The equivalent circuit of a NSN transistor, 
(b) The density of state when there are odd number of 
electrons in the superconducting island. 

Figure 1 (a) shows the equivalent circuit of 
the NSN transistor. A superconducting island ex- 
changes quasiparticles (QPs) with a left (right) 
lead via a tunnel junction characterized by the 
tunneling matrix element 7l(r) and is coupled to 
a left (right) lead and a gate via capacitors Cl(r) 
and Cq. In the following, we consider the symmet- 
ric system where Tl = Tr and Cl = Cr. We limit 
ourselves to the discussion on charge fluctuation 
at zero temperature where the lowest two energy 
states make a dominant contribution. In the lan- 
guage of effective spin cr, the total Hamiltonian is 
written as 

-H" = ^ Srfe alkn^rkn + ^ (T^e*'^''^ a|j.„arfe'«o-+ 



h.c. 



Ao(o-2 



r=L,R 

k,k' ,n 

l)/2 + So, 



(1) 



where the first term is the Hamiltonian of QPs in 
the lead r (=L,R) and the island (r=I). a,-kn is the 
annihilation operator of a QP with wave vector 
k, transverse channel n which includes the spin. 
We consider the limit of large number of transfer 
channels, A^ch- We assume that an unpaired elec- 
tron occupies the lowest energy level above the su- 
perconducting gap. Figure 1(b) shows the normal- 
ized density of state (DOS) in the island pi{e)/pi. 
The second term describes the tunneling across 
junctions and simultaneous change of the charge 
state, if is the phase difference between the left 
and right leads and kl = — kr = 1/2. The third 
term describes the excitation energy, namely the 
energy difference between two charge states, which 



is written with a gate charge Qq as Ag = Ec(l — 
2QG/e) + A. The last term Eq = EciQc/e)^ is 
the charging energy for the charge state |0) where 
there are even number of electrons in the island. 

We employ a mapping of the spin-1/2 operator 
onto two fermion operators, c and d, [9] and use the 
Schwinger-Keldysh approach [8, 10]. The approxi- 
mate generating functional is expressed by T4^ = 
— zftlnTr[C~^], where trace is performed over the 
time t defined on the closed time-path C which 
goes from — oo to oo, returns along the real axis 
and closes at — oo — ih/3. Here Gc is the full Green 
function of c field defined as G~^(ti, ^2) — (ifidt — 
h(ti))d{ti, ^2) — Sc(ti, ^2), where S is (5-function de- 
fined on C and we replace Aq with h{t) in order to 
preform the functional derivative. h(t) is regarded 
as the "external magnetic field" acting on the spin 
a. The c field self-energy is given by the product of 
particle hole propagator and the free Green func- 
tion of d field: Ttc(ti,t2) = ihg^(tl,t2)Q.{tl,t2)■ 
The average charge in the island is calculated by 
functional derivative of the generating functional 
W with respect to the magnetic field, Q/e — 1/2 ~ 



, where /ia(c) is the rela- 



<5W^/<5/lA(i)L,=0...=Ao 

tive (center-of-mass) coordinate[10] of the mag- 
netic field. The phase difference is fixed as (fdt) = 
eVt/h and ip/^ (t) =0. The average charge is rewrit- 
ten as. 



Q/e- 



de a^(£) 



ImGf(£), 



(2) 



where a^ (e) = -~2'Ki X]r=L r <^r Ip(^ — A — K^eV) \ 
and a^{e) ~ ^'''^ Sr=L r "^r P(^ ^ A — K,-eV) 
are the Keldysh and retarded component of the 
particle-hole propagator, respectively. Here a° = 
NchT'^PiPr and p^ is DOS in the lead r. The spec- 
tral density of p(e) is suppressed around e = due 
to the superconducting gap [8]. 

In Eq. (2), the imaginary part of the retarded full 
c field Green function Cf (e) ^ l/(e - Aq - Sf (e)) 
represents the spectral function of the excitation 
of the charge state[8]. Im[S]^] — —7 is the life-time 
broadening caused by dissipative charge fluctua- 
tion. The real part of the self-energy is the shift of 
a pole of the Green function. It is noticed that at 



zero-bias voltage Eq. (2) is formally equivalent to 
the expression obtained within the resonant tun- 
neling approximation for normal metal island [2]. 
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Fig. 2. (a) The average charge as a function of the excita- 
tion energy at zero-bias vohage for qq =0.1 (solid line), 0.2 
(dashed line) 0.3 (dot-dashed line) and 0.4 (dotted line). 
Vertical lines indicate positions of degeneracy points [6]. 
Two insets show schematic diagrams of dominant process 
of charge fluctuation in CB regime for odd state (panel 
(a-1)) and that for even state (panel (a-2)). In the panel 
(a-l)/(a-2), the island chemical potential is under/above 
that of leads. In CB regime for odd state, the charge fluc- 
tuation is blocked due to the superconducting gap. (b) The 
ground-state energy as a function of the excitation energy 
for various ao ■ The long dashed lines indicate the charging 
energy of even (Eq) and odd (-Ei) state, respectively. As qq 
increases, the local maximum point shifts leftwards. The 
inset of the panel (b) shows a schematic diagram of gate 
charge dependence of the charging energy. Solid and dashed 
lines indicate the energy with and without charge fluctua- 
tion, respectively. The energy correction in CB regime for 
even state is larger than that for odd state. Parameters 
are chosen as A/Ec = 0.1 and p^Ec = 10^^. The value 
ApJ* = 10^ is the same order as an Al island whose volume 
is lO'^nm^. 

Figure 2(a) shows the average charge as a func- 



tion of the excitation energy at zero-bias volt- 
age for various dimensionless conductance oq = 
X]r=L R "^r ■ Vertical Hues indicate positions of 
degeneracy points. Each position is given by the 
real part of the pole of the full c field Green func- 
tion. We can see the degeneracy point is shifted 
leftwards as ao increases. The deviation of the 
average charge from the discretized value, in CB 
regime for even state is larger than that for odd 
state. This asymmetric behavior around the de- 
generacy point can be understood by considering 
processes of charge fluctuation. Two insets show 
the schematic diagrams of the dominant process of 
quantum fluctuation of QP in CB regime for odd 
(Fig. 2(a-l)) and even (Fig. 2(a-2)) states. In the 
former regime, the dominant process corresponds 
to a particle-hole excitation with a particle in one 
lead and a hole under the superconducting gap in 
the island (Fig. 2(a-l)). In this case, a low-energy 
particle-hole excitation, which makes a large con- 
tribution to the charge fluctuation, is blocked due 
to the superconducting gap. On the other hand, 
in CB regime for even state, the dominant process 
corresponds to a particle-hole excitation with a 
hole in one lead and particle in the island (Fig. 
2(a-2)). Because a low-energy particle-hole excita- 
tion is enhanced due to the large DOS at the edge 
of the superconducting gap, the charge fluctuation 
is large. 

The excitation energy dependence of the ground- 
state energy for various ao is plotted in Fig. 2(b). 
The ground-state energy is expressed as 



E — En 



d£-Imln(Gf(e)) 



(3) 



The energy correction in CB regime for even state 
is larger than that for odd state, because the charge 
fluctuation is large in the former regime. As ao in- 
creases, the local maximum point shifts leftwards 
and the degeneracy point shifts in the same direc- 
tion. This behavior can be explained by consider- 
ing the charging energy renormalization[8]. The in- 
set of Fig. 2(b) shows a schematic diagram of gate 
charge dependence of the charging energy. Solid 



and dashed lines indicate the charging energy with 
and without charge fluctuation, respectively. One 
can see that the energy correction caused by the 
charging energy renormalization in CB regime for 
even state is larger than that for odd state. 
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Fig. 3. The 3D plot of the average charge versus the 
excitation energy and the bias voltage for (a) qq =0.1 
and (b) 0.3. The kink structure is seen at the boundary 
of CB regime for even state. The schematic diagrams of 
dissipative charge fluctuation process at the boundary of 
CB regime for odd state (c) and that for even state (d). 
In the panel (c)/(d), the island chemical potential matches 
that of the left/right lead, (c) The fluctuation is blocked 
by the superconducting gap. (d) The fluctuation increases 
with increasing bias voltage. 

The average charge at finite bias voltage is shown 
in Figs. 3(a) and 3(b). Even for large ao (Fig. 
3 (b) ) , there is kink structure at the boundary of CB 
regime for even state. The kink structure is robust 
against charge fluctuation in the range of e^ < 2A, 
because dissipative charge fluctuation is blocked 
due to the superconducting gap (Fig. 3 (c)), and 
only the unpaired electron contributes to the life- 
time broadening: 7 ~ 7rao/2p°. On the other hand 
at the boundary of CB regime for odd state, the 
dissipative charge fluctuation increases with in- 
creasing bias voltage (Fig. 3 (d)), and the life-time 
broadening 7 is proportional to the number of QP 
states between the chemical potential of the left 
lead and that of the island: 7 ^ TraoVA eV. There- 
fore, a small variation of the average charge at the 
boundary is smeared by applied bias voltage. 

In conclusion, we have theoretically investigated 
the charge fluctuation between even and odd states 



of the superconducting island. The charge fluctua- 
tion is suppressed due to the superconducting gap 
in CB regime for odd state, while it is enhanced 
in CB regime for even state due to the large DOS 
at the superconducting gap edge. The large charge 
fluctuation causes the large energy correction in 
CB regime for even state, which can be under- 
stood by considering the charging energy renor- 
malization. At finite bias voltage, there is the kink 
structure in the 3D plot of the average charge at 
the boundary of CB regime for even state, even 
for large ao . This is because the dissipative charge 
fluctuation is suppressed due to the superconduct- 
ing gap. 
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